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A potential implementation of quantum-information schemes in semiconductor nanostructures is 
studied. To this end, the formal theory of quantum encoding for avoiding errors is recalled and the 
existence of noiseless states for model systems is discussed. Based on this theoretical framework, we 
analyze the possibility of designing noiseless quantum codes in realistic semiconductor structures. In 
the specific implementation considered, information is encoded in the lowest energy sector of charge 
excitations of a linear array of quantum dots. The decoherence channel considered is electron-phonon 
coupling We show that besides the well-known phonon bottleneck, reducing single-qubit decoherence, 
suitable many-qubit initial preparation as well as register design may enhance the decoherence time 
by several orders of magnitude. This behaviour stems from the eífective one-dimensional character 
of the phononic environment in the relevant region of physical parameters. 
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I. INTRODUCTION 



Devices using unique quantum-mechanical features can 
perform information processing in a much more efïicient 
— or even unattainable — way than those relying just on 
classical physics. This fundamental discovery has stimu- 
lated in the last few years a big deal of work and scientific. 
debatés in the new born field of Quantum Computation.EI 
From a conceptual point of view these results represent 
a serious challenge to the time-honored notion of univer- 
sal computational schemes independent of an underlying 
physical theory: information as well as computation are 
intrinsically physical. On the other hand, physical real- 
izations of a Quantum Computer would result in tremen- 
dous practical advantages. 

The key ingredients which endow QC devices with 
computational capabilities that supersede their classical 
counterparts are basically: (i) the linear structure of their 
state space; (ii) the unitary character of their dynam- 
ical evolution; (iii) the tensorised form of multiparticle 
state spaces. The first two properties allow for a paral·lel 
processing of an arbitrary number of data sets, encoded 
in suitable quantum states. By resorting to quantum 
interference, between different computational branches, 
one can selectively amplify desired parts of the state vec- 
tor in order to optimize the probability that a final (i.e., 
read-out) measurement will give us the information we 
were looking for. Point (iii) represents another striking 
departure from classicality: due to entanglement, com- 
bining different quantum systems results in an exponen- 
tial growth of the availablc coding space; moreover, the 
tensor-product structure is at the very basis of many ef- 
ficient quantum manipulations. 



Unfortunately, all this holds just for closed quantum 
systems. Real systems are unavoidably coupled with en- 
vironmental (i.e., non computational) degrees of freedom. 
Such open character spoils points (i) and (ii) eventually 
turning quantum computing to classical. Different com- 
putational branches get entangled with different (quasi- 
orthogonal) quantum states of the environment and their 
interference is then no longer observable. From a math- 
ematical point of view, the relevant state space, given 
by density matrices, has now a convex structure and the 
allowed quantum dynamics is described by CP-maps.EI 
Initial pure preparations are typically corrupted on ex- 
tremely short time-scales due to quantum-coherence loss 
that makes them mixed: the initial information irre- 
versibly leaks out from the system into the huge number 
of uncontrollable degrees of freedom of the environment. 
This nhenomenon — the so called decoherence problem 
in QCli — represents the major obstacle for the exper- 
imental realization of any quantum-computing system. 
Other challenging requirements are of course given by 
the necessity of being able to perform on a system, with 
a well-defined state space, long coherent quantum ma- 
nipulations (gating) , precise quantum-state synthesis and 
detection as well. 

A major theoretical achievement has been made by 
showing that one can, in principle, actively stabilizc 
quantum states by means of Quantum Error Correction.Q 
The latter, built in analogy with its classical counterpart, 
assumes that the quantum bits (qubits) are coupled to 
independent environments. The information is then en- 
coded in a subtle redundant way that allows, monitoring 
the systems and conditionally carrying on suitable quan- 
tum operations, to tolerate a certain (small) amount of 
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decoherence and imperfect gating as wellü 

It is basically the need of dealing with systems suf- 
ficiently decoupled from the external environment that, 
up-to-now, has limited the existing realizations to atòmic 
and molecular implementations. Furthermore, the ex- 
tremely advanced technological state-of-the-art in these 
fields [-allows for the manipulations required in simple 
QC'sü However, any interesting QC would require a large 
number of quantum gates and qubits as well, and all the 
present approaches suffer from the problem of scalability 
to large, i.e. highly integrated, quantum processors. 

One is then naturally led to consider the viability of 
solid-state implementations. In particular, by resorting 
to present semiconductor technology, one might bene- 
fit synergeticaliy n from the recent progress in ultrafast 
optoelectronicaïíLand in nanostructure fabrication and 
characterization.E2l 

The first drawback of such a kind of proposal is that 
the typical decoherence time Trj in semiconductors is of 
the order of picoseconds. On the other hand, the relevant 
parameter is the ratio between the typical time-scale of 
gating tq and tq. Roughly speaking, td/tq represents 
the number of elementary (coherent) operations that one 
could perform on the system before its coherence being 
lost. j_. 

DiVincenzo and LossEil have proposed to use non- 
cquilibrium spin dynamics in quantum dots for quantum 
computation. This exploits the low decoherence of spin 
degrees of freedom in comparison to the one of charge 
excitations, being the former much less coupled with the 
environment. Nevertheless, the required magnètic gat- 
ing is extremely challenging from a technological point of 
view, and the ratio td /tg does not allow for the number 
of gate operations within the decoherence time required 
by concrete QC's. 

Ultrafast làser technology is now able to generate elec- 
trònic excitations on a sub-picosecond time-scale and 
to perform on such states a variety of coherent-carrier- 
control operations .0 If one can speculate to resort to such 
a technology for realizing gating of charge degrees of free- 
dom then coherence times on nano/microsecond seales 
can be regarded as "long" ones. 

In this paper we analyse in a detailed way the re- 
cent idea_of implementing Quantum Error Avoiding 
strategies.Ej The goal here is to suppress decoherence in 
a quantum register realized by the lowest energy charge 
excitations of a semiconductor quantum-dot array.1121 In 
this case, the noise source is given by electron-phonon 
scattering, which is recognized to be-Üie most efficicnt 
decoherence channel in such a systenocl. 

Despitc of the a priori complexity of the three- 
dimensional (3D) phononic environment, we will show 
that the underlying dynamical-symmetry allows, by 
means of a proper quantum encoding, to increase the 
decoherence time by several orders of magnitude with 
respect to the bulk value. The focus of the present paper 
is mostly conceptual and the problem of actual prepara- 
tion/manipulation of the resulting codewords will not be 



addressed. 

The paper is organized as follows. In Sect. |ïï] the for- 
mal theory of subdecoherent quantum encoding is pre- 
sented and discussed. Section III deals with the appli- 



cation of the proposed subdecoherence theory to real- 
istic semiconductor-based nanostructures; More specif- 
ically, we will choose as quantum register an array of 
semiconductor quantum dots and for this particular sys- 
tem we will study the potential sources of decoherence. 
In Sect. IV we shall present a detailed investigation of 
decoherence in our quantum-dot array. In addition to a 
short-time analysis, we will present time-dependent simu- 
lations corresponding to a numerical solution of the Mas- 
ter equation. They will show that by means of a proper 
initial many-electron state preparation it is possible to 
extend the carrier-phonon decoherence time up to the /is 
seale. Finally, in Sect. ^ we will summarise and draw 
some conclusions. Appendix A is devoted to a formal 
analysis of the so-called Circular Model, which will turn 
out to play a major role in the semiconductor-based im- 
plcmcntation considered. 



II. THEORY OF SUBDECOHERENT QUANTUM 
ENCODING 



In this section—we recali the bàsics of the theory of 
Noiseless Coding}3 in the framework of a Master Equa^ 
tion (ME) formalism, for the register subdynamics.o 
Generally speaking, these strategies for preserving quan- 
tum coherence rely on the possibility to design an open 
quantum system 1Z in such a way that i) the environment 
£ is effectively coupled only with a subset of the degrees 
of freedom of 1Z. Information is then encoded in the por- 
tion C of Hilbert space spanned by the remaining (decou- 
pled) degrees of freedom, ii) The environment is coupled 
to subset of states C in a state independent fashion. In 
both cases £ is not able to extract information from C : 
the quantum coherence is then passively stabilized. From 
the above points it should be clear the first and major 
departure from the Error Correction paradigm: here one 
assumes the environmental noise to be correlated. £ is 
coupled, in a strongly state- dependent way, with collective 
states of S. 

Before embarking in a detailed analysis of sub- 
decoherence let us shortly discuss two very simple ex- 
amples, that show how this notion can come about. 

i) Let us consider TV isospectral linear oscillators H-r = 
uj ^2j—i bj bj coupled with the vacuum fluctuations i.e., 
zero temperature, of a bosonic field afc by an Hamilto- 
nian of the form Hj = X^7c(fffcj a Í bj +h.c). Suppose now 
that gkj — gk Vj. By introducing the Fourier transformed 
operators b q = l/V'TVX^ & %q ^bj (bosons as well) one im- 
mediately sees that only the zero-modes are actually cou- 
pled: Hi = &o(E*0* a *) + h-C- and H n = vJ2q b q b q- 

Therefore, any state of the (infinite-dimensional) sub- 
space 
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C = |0) o (g)W, C 1 ) 

<?>o 

will evolve unafïected by the environment in that Hj C ® 
|0) £ =0. 

ii) Let the system-environment interaction Hamilto- 
nian be of the form Hx = ^2 ® E^, where e 
End7Y i ^(A' = 1Z, £), Moreover, let us suppose that 
the Hermitian R^s are commuting operators, i.e., they 
span an abelian àlgebra A. Let C C a simultaneous 
eigenspace of .4. This means that 

#r|c=X>p®£„ (p^eR), (2) 

in other words, if one restricts himself to C the interac- 
tion with the environment amounts simply to a state- 
independent renormalization of Hg. It is then clear 
that — provided C is invariant under the system self- 
Hamiltonian H-r — any initial preparation in C evolves in 
a unitary fashion regardless the strength of the system- 
environment coupling and the environment initial state 
as well. Of course, for all this to be usefïïl in quantum 
encoding one must have dimC > 1. 

A. Master-equation approach 

The system under investigation TZ is given by N iden- 
ti cal two-level systems (N-qubits quantum register ), rep- 
resenting our computational degrees of freedom, couplcd 
with an external (uncontrollable) environment. The reg- 
ister TZ will be described in the spin 1/2 language by 
means of the usual Pauli spin matrices {of o\ }{Li gen- 
erating N local sl(2) algebras 

[af,aJ}=2 5 ij o-i : [of , of] = ±%o±. (3) 

The collective spin operators S a — J2iLi a ? > ( a = z ) 
span a sl{2) àlgebra as well, it will be referred to as the 
global sl(2). The environment £ will be described by a 
set of non interacting harmònic oscillators with bosonic 
field operators [b' k , = Skk'- 

The total Hamiltonian is assumed to be H = H-o + 
H e + Hi, where H n = E S z and H £ = £ fc w fc &Iò fc 
are, respectively, the register and the environment self- 
Hamiltonians. Here, E represents the energy spacing 
between states \0)i and \l)i in cach qubit. The TZ — £ 
interaction is given by 

£Fr = X>w&Ur + h. C ). ( 4 ) 

ki 

Let us now briefly recali the Standard Born-Markov 
scheme for tracing out the £ degrees of freedom and 
obtaining a Master equation for the register subdynam- 
ics. The Liouville-von Neumann equation for the total 
density matrix of TZ ® £ in the interaction picture reads 



i dtp = [Hx, p\. One assumes a factorized initial condi- 
tion j5(0) = p ® 0. After a formal time integration one 
obtains 

p(t)=p(0)+ f oIt[Hx{t),p(t)] 
Jo 

= p(0)~i í dr[Hx(r),p(0)} 
Jo 

+ {-if fdr i dr'[Hx(r),[Hx(r'),p(r')}} (5) 
Jo Jo 

Now we set p(r') = p{r') ® íl (Í2 ~ e~P He ) and we per- 
form a partial trace over £ in order to get an equation 
for the reduced density matrix of TZ: p(t) = tr f p(t). The 
resulting ME is of the form p = C(p). The Liouvillian 
superoperator C is given by the sum of two contributions: 
C u representing the unitary component of the dynamics 
ruled by the (renormalized) register self-Hamiltonian; Ca 
describing the irreversible decoherence/dissipation pro- 
cesses induced by the coupling with the external bath. 
By denoting with H-r the unperturbed register self- 
Hamiltonian, one has that C u (p) — i/Ti[p, H-r + SHn] 
where the environmcnt-induccd SH-r is given by 

r)=± ij = l 

These contributions — usually referred to as the Lamb- 
shift terms — describe a sort of qubit-qubit effective in- 
teraction mediated by the external environment. The 
dissipative Liouvillian is given by Cd = X)r;=± ^2' wnere 

1 N 

c » = 2ï E r S (w p> -7 r 'i + [°i> p°J n ]) . (7) 

ij=l 

Here, the term rj — — (77 = +) is associated to deex- 
citation (excitation) processes of the qubits by emission 
(absorption) of bosonic quanta. The Hermitian matrices 
r and A are the input data defining our ME, their ac- 
tual form depends on the details of the physical constants 
{E, {u>k}k {dki}, etc.) and will be given later. 

As far as the analysis of this section is concerned is suf- 
ficient to know that T > 0. One can go on with general 
considerations by diagonalizing in order to obtain 
the jGanonical form for the dissipative part of the Liouvil- 
lianlia 

1 N 

Cd ^=2h E X M\ L >> L--"] + [Ll,pL-*]), 
»j=±,fí=i 

(8) 

where {A^} are the (non-negative) eigenvalues of r^. 
Moreover, L 1 ^ = X)i M f a l.i u i denoting the components 
of the eigenvectors of T^. The L^'s will be referred to 
as the Lindblad operators. The operator (Lie) àlgebra A 
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spanned by the Lindblad operators contains the informa- 
tion about the existence of coding spaces stable at least 
on a short-time-scale. The finite-time stability depends 
on the interplay between the dissipative and the unitary 
components of the Liouvillian in Un . 

In order to quantify the efficiency of the environment 
in destroying quantum coherence it is useful to define a 
(first-order) decoherence time (rate) t\ (t^-,) by means 
of the short-time expansion of the fidelityíÉ for pure ini- 
tial state preparations 

F(í)^(V|p(í)|V) = l-^+o(í 2 ). (9) 

From Eq. (@) one obtains 

N 

rï l M]= E K(\\Lim 2 -ML*\i;)\ 2 ). (10) 

This expression is nothing but a sort of fluctuation- 
dissipation relation connecting the dispersion of the Lind- 
blad operators U 1 ^ in the initial register state with the 
rate at which quantum coherence is destroyed. It is im- 
portant to point out that the unitary component of the 
Liouvillian does not contribute to the first-order deco- 
herence time. If t 1 ~ 1 [|'í/')] = then the state \ip) will be 
called subdecoherent and a linear subspace C C Hn will 
be referred to as a subdecoherent code. 

In general, the register Hilbert space splits in A- 
invariant subspaces, 

■Hn = @j® n r i l U { ; , \ (11) 

where J labels the irreducible representations (irrep) of 
A, and the integers nj are the associated multiplicitics 
[Tír^ = Tí),, .] The singlet sector C of A is the direct 
sum_Xpossibly empty) of the one-dimensional irreps. In 
Ref.EJ it has been shown that, for non-abelian A, the 
sub-decoherent codes coincide with C. In an equivalent 
group-theoretic language one can say that the code C 
is the subspace of vectors invariant under the action 
of group G — exp A generated (infinitcsimally) by the 
Lindblad operators L^. [C is the trivial (J-representation 
space.] This group acts, of course, on the general mixed 
states: p^l pX >, (X e Q). The same argument holds 
for the subdecoherent (pure) states. When C is invari- 
ant under the action of H^; then the contribution to the 
dynamics of Cd vanishes 

p= \i>){ií)\ h-» e - itH '«-pe iiH * (V|V) G C, t > 0) (12) 

The finite-time evolution is unitary, in this case C will 
be referred to as noiseless code: quantum coherence is 
preserved - in principle - for an arbitrarily long time. 
When C is not ií^-invariant the initial preparation 
on a greater time-scale leaks out from the code and its 
quantum coherence will be eventually washed out. For 
instance, the condition [H' n , A] = 0, sufhces to have such 



a noiseless coding or even that H' n belonging to the as- 
sociative operator àlgebra A a generated by the L^s and 
the identity operator. Notice that if C is subdecoherent 
for the L M 's it is subdecoherent for any set of Lindblad 
operators included in A a - 

From a physical point of view, the àlgebra A of Lind- 
blad operators represents the set of the register modes 
that are incoherently excited by the environment; look- 
ing for states that are annihilatcd by as many Linbdlad 
operators as possible is then as looking for states that 
are " vacua " for the largest number of such excitatiopa 
and therefore maximally decoupled with environment x3 
It is important to emphasize that such a decoupling can 
be achieved thanks the algebraic-dynamical structure of 
the model without any assumptions about the (weakncss) 
of the register-environment interaction. Loosely speak- 
ing, one can say that for genèric T's, the Liouvillian is 
such that, given any register preparation, the environ- 
ment forces the coding system to explore the totality of 
its Hilbert space so that there is no safe place where stor- 
ing quantum information, instead for some "màgic" T the 
Lindblad àlgebra gets smaller allowing just for a limitcd 
probing of the register space of states by the environment 
strongly dependent on the initial register data: free room 
is left for " hiding " quantum information. 

Rather interestingly, the problem of analysing state 
stability against decoherence can be cast in a Hamilto- 
nian form by observing that, for an initial condition \tp) 
that is a S z -eigenstate one has t^ 1 = {ip \ H where 

N N 

È = E w^ = E( r ^+f^) (13) 

In other words: the problem of finding decoherence rates 
is mapped onto the spectral problem for the (positive) op- 
erator (jI3[). In particular, "robust" states (i.e., the ones 
with minimal decoherence rates) are ground states of H. 
Let En denote the lowest eigenvalue of H. En = 
means that there exist sub-decoherent states, in this 
case C = Ker H and dN = dim Ker H gives the dimen- 
sion of the code. The sub-decoherence property is sta- 
ble against small perturbations of the state. Indeed if 
G Kerií h-> + \6tp) then írf 1 = {Stp\ H \5%j)) > 



B. A Simple Example 

To better illustrate the situation let us consider the 
N = 2 case. The model ([L3h is soluble in elemen- 

tary way0. We assume r£f = = and 

r if ) = r 2i ] = r(±) A moreover > rW . From posi- 
tivity it follows that \f3\ < 1. The spectrum is given by 

£ t)S = (r^ + r(+))(i±/3), (14) 
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with eigenstates given respectively by 

|11), |00>, 2-Vafloi) ± |10)). 

if r(+) > o, for \p\ < (r(-) - r(+))(r(-) + rM)- 1 = p c 

one has E 2 = Eq , for (3 > f3 c (/? < (3 C one finds E = E s 
(E = E t .) r(+) = =>■ ^2 = £ o = 0. Finally for 
(3 = ±1 one has again i?2 = 0, with eigenstates given 
by IV's.í)- In summary, subdecoherent states exist in a 
subsct of the boundary of the T manifold. This result 
is quite general: for genèric T's one has £W[r] > 0, the 
subdecoherence condition -EjvfT] = is fulfilled just in a 
"zero-measure" set of the Hamiltonian models (^3|). Of 
course this is simply due to the fact that for a genèric 
T's gives rise to a Lindblad àlgebra A is too large for 
admitting a (non-trivial) singlet sector. 

Turning back to the general N case, to exemplify the 
collectivc nature of the decoherence-dissipation dynamics 
let us consider the states (TV even) 

\^ ym ) = (sY /2 \o), 

|Z>7> = ® W )6i>(|01) - (_i)-y(*.i)|io» ii , (15) 

where D is a dimer partition of the qubit array, and 
7: T> — ► {0, 1}. The first state in ( |Ï5| ) is simply the totally 
symmetric S z — state (belonging to the sl(2) multiplct 
of the vacuum) whereas the \T>j)'s are products of sin- 
glet or triplet pair-states depending on the signature 7 
of the register dimer partition T>. This latter family of 
states (Q~5Ï) will play an important role in the following. 
Notice that, for 7 = 0, one gets global sl(2) singlets cor- 
responding to zero total angular momentum S 2 . In terms 
of Hadamard transformations and controlled-not opera- 
tors the |2?<r)'s can be sinthetized as follows from a pure 
product state 



\Vj) = (8)( /jm ) eP cnot / „ l ií ; m) + 1, 1] 



Ini 



(16) 



With a straightforward calculations one finds that the 
first-order decoherence rates of states (|1]) are given re- 
spectively by (tq/to) -1 = f a (T), (a = sym,Vj) in 
which Tq 1 ~Tq N/2 is the decoherence rate for uncorre- 
lated qubits and (Tíj = Tij/To) 



fs: 



N 



/ ü7 = i--k j2 (-iy {i ' j) ïij 



(17) 



where To = Tu, (i = 1, ... , N) The / Q 's contain the in- 
formation about the degree of many-qubit correlation in 
the decay process: if T oc I one has f a = 1 the qubits 
decohere independently. 



III. APPLICATION TO SEMICONDUCTOR 
NANOSTRUCTURES 



to realistic, i.e., state-of-the-art, semiconductor-based 
nanostructures. Since in semiconductors the primary 
source of decoherence is known to be carrier-phonon scat- 
tcring, we will consider as prototypical systems quasi 
zero-dimensional (0D) structures, for which the reduccd 
phase-space available allows for a significant suppression 
of phonon-induced energy relaxation and dephasing. 

We will choose as prototype of quantum register an 
array of semiconductor quantum dots. In particular, 
we will consider as quantum dot (QD) a GaAs/AlGaAs 
structure similar to that studied in [Ï8|. Here, vari- 
ous effects due to carrier-carrier interaction will not be 
considered. This is, of course, a potential limitation 
of our analysis, especially in relation to state prepara- 
tion/manipulation (not addressed in this paper). Indeed, 
the latter requires a controllable source of entanglement, 
i.e., a qubit-qubit interaction thatjpight be provided by 
"switchable" Coulomb couplings.Ej On the other hand, 
our coding states will involve single-electron occupations 
only; For such states the intra-dot Coulomb repulsion 
is clearly absent, while the inter-dot one at the dis- 
tances relesant for our quantum encoding is found to be 
negligible .cJ Moreover, since the system under considera- 
tion is based on intrinsic III- V materials, carrier-impurity 
scattering is negligible. 

Generally speaking, Hamiltonian modifications will re- 
sult in leakage from the coding subspace only on a longer 
time-scale with respect to the phonon-scattering one, i.e., 
it does not affect the stability classification based on 
t\ (see Sect. |J). Finally, we would like to stress that 
there existe a whole class of interactions leaving the code 
invariantJ13 



A. Free-Carrier States in the Quantum-Dot Array 



The confincment potential V 0I) giving rise to the quasi- 
0D carrier states in such a QD structure is properly 
described in terms of a quantum-well (QW) profile V" 
along the growth direction of the structure plus a two- 
dimensional (2D) parabòlic potential V in the nor- 
mal plane. More specifically, a carrier within the i-th 
QD structure is described by the following single-particle 
Hamiltonian 



2m* 



+ V m (r)= I - 



2m* 



2m* 



+ VÍ(r 



iíll 



(18) 



where 



In this section we shall discuss a potential application 
of the above sub-decoherent quantum-encoding strategies 



^(r^H-mVlrJ 2 , 



(19) 
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is the 2D harmonic-oscillator potential in the (x, y) plane 
perpendicular to the (z) array axis (which coincides with 

the growth axis of the QD structure), while V i (r||) is 
a 1D square-well potential centercd at rj| = iai with 

width d and infinite walls0 a being the array period- 
icity, i.e., the inter-dot distance. This choice for the 
single-particle Hamiltonian, even though not generally 
vàlid, well describes the OD carrier confinement of the 
low-energy states in the QD structure, which are the only 
relevant states for the quantum encoding considered. We 
would like to point out that the very same QD model 
turned out to be able to explain, in-a quantitative way, 
the addition spectra reported in iïq.E3 

The Hamiltonian ([18] )is elementary soluble, its spec- 
trum being the sum of the parallel and perpendicular 
contributions: 

= EÍ + EI = {n x +n y + l )fwJ + — — -r . 20 

2m d z 

The corresponding 3D eigenstates will be factorized ac- 
cording to: 

<t>i,nv{v) = <t>t, ny {vs_)(t>Ur\\ -ia). (21) 

The total free-carrier Hamiltonian describing our QD ar- 
ray can then be expressed in the (second- quantized) form 

Hn = ^ e " c L c "*> ( 22 ) 

i,ot 

where the fermionic operators c\ a ( Cj a ) create (destroy) 
an electron in the z-th QD in state a = n x n y v. 



C. The Qubit Register 

In the proposed information-encoding scheme the sin- 
gle qubit is given by the two lowest energy levels of the 
QD structure. Since the width d of the GaAs QW region 
is typically of the order of few nanometers, the energy 
splitting due to the quantization along the growth di- 
rcction is much larger than the confinement energy tko 
induced by the 2D parabòlic potential V (typically of 
a few meV). Thus, the two lowest-energy states — state 
|0) and |1) — realizing our qubit are given by products of 
the QW ground state times the gmund or first excited 
state of the 2D parabòlic potentiaÈJ. 

More specifically, they are given by 

(r\l) l = ^(x) ( f > i(y)4 (z) (26) 

where 

- C e- a " x \ C = (2a /7T) 1 / 4 , a = ^ 
4>Í (s) = Ci x e- ao x ' , Ci = 2 al /4 (2/tt) I ' i (27) 

are, respectively, the ground and first excited states of 
the harmònic oscillator in the (perpendicular) xy plane, 
and 

4 >0 (z) = C z cos[^(z-ia)], ,C Z = V2/d (28) 



B. Carrier-Phonon Coupling 

The Hamiltonian describing-the free phonons of a semi- 
conductor crystal is given byCJ 



Aq 



(23) 



where A and q denote, respectively, the phonon mode 
(e.g. acoustic, optical, etc) and the phonon wavevector. 

The coupling of phonons with the electrons in the QD 
array is described by the following carrier-phonon inter- 
action Hamiltonian: 



H I= [9i a ,i'cx';Xq_ c L· b Xq C i'c t ' + h - c - 

ia,i'a';Aq 

Where 



9i 



Aq = 3A q / ^ Q (r)e lqr ^ Q ,(r)dr 



(24) 



(25) 



are the matrix elements of the phonon potential between 
the quasi-OD states ia and i'a! . The explicit form of the 
coupling constant gx^ depends on the particular phonon 
mode. 



is the ground state of the i-th quantum-well potential 
parallel to the array axis {4>\ (z) — for \z — i a\ > d/2). 

Notice that the only dependence on the QD label i of 
the qubit states is in the z-component of the wavefunc- 
tion. 

Since we are restricting ourselves to the low-energy sec- 
tor a = 0, 1 in the absence of inter-dot (i ^ i') transi- 
tions, the only relevant fermionic bilinears in Eq. (|24|) are 
given by Xi — c^Cjo and their conjugates. Consistently 
with the commutation relations [Xi, xj] = Sijfaj—n®) = 
2af, these bilinears can be described by the spin 1/2 op- 
erators of. Let |0) = nï=i c !o l vac ) the reference state 
built over the electron vacuum by occupying all the |0)j. 
Our reduced Hilbert space containing the computational 
degrees of frecdom is then given by 



iY 



JV 



H* = span{JjA-f |0) | Oí = 0, 1} 



(29) 



i=l 



Any process inducing transitions out of this subspace will 
result in a computational error. Let A being the energy 
gap between |1) and the highcr excited states (in the 



G 



present case A = fiu) and T the environment (i.e., lat- 
tice) temperature; this sort of leakage errors occur with 
low probability as long as A ^> ks T. 

By denoting with E = e^i — e^o — the energy spac- 
ing between our two qubit levels, the free-carrier Hamil- 
tonian ( ^2| ) for our qubit register, i.e. restricted to the 
low-energy sector a — 0, 1, can then be written as 



JV 



where cr| denotes the usual diagonal Pauli matrix acting 
on the i-th qubit. 

Let us now consider again the carrier-phonon interac- 
tion Hamiltonian (|2^). Within the carrier model con- 
sidered, wavefunctions corresponding to difïerent QD's 
do not overlap; thus one has gia.i'a'-.Xq = for i ^ i', 
i.e., phonons induce intra-dot (intra-qubit) transitions 
only. The coupling constants associated to the relevant 
elementary processes in our qubit register are <7i,Aq = 
5ïi,io ; Aq,5i,Aq = 3io,ii;Aq- More specifically, starting from 
the explicit form of the single-particle wavefunctions 
<f> in (||) one finds g iyX q = 9Xq9x(qx) 9y{ly) 9z[lz, *i) 
[q = ÍQx, q y , qz)] where 



9x{qx) 
9iz(qz) 



)=exp(--í^) 



sa 



2 a, 



1 t q v ^ 



8tt 2 sin(a z d/2) 
d 3 q z ql - q 2 z 



z 



* (31) 



and qo = 2n/d 

Within these assumptions the carrier-phonon interac- 
tion Hamiltonian ( ^4|) can be cast in to the form (Ü|): 



Hi = ^{g^blal 



+ h. c) , 



(32) 



ki 



where the bosonic label k now corresponds to the phonon 
modes of the crystal, i.e., k = Aq. 

Following the Born-Markov procedure discussed in 
Sect. 0, one finds the Jbllowing result for the matrices 
r, A, defining our MEEl 



V ij = 2 n 9^ 9kj ( n k + 9(t) ) S(hiü k - E), 



9ki 9k 3 

huj k - E 



(n k + 9(T)) 



(33) 



Here, 9 is the customary Heaviside function, and V de- 
notes the principal part. From these relations it follows 
that and A^ are hermitian as expected. Further- 
more > and T^' > rW. Since for the QD struc- 
tures considered the energy splitting E is typically much 
smaller than the optical-phonon energy [36 mev in GaAs] 
the only phonon modes k = Aq involved are the acoustic 



ones. In this case, by considering carrier-phonon coupling 
due to deformation potcntial, one has g(q) = \J 2pVc ' 
where e is the scalar lattice deformation, p and V the 
crystal mass-density and volume, while c is the sound 
velocity. 

Let us now focus on the explicit form of the function 
T in íil i.e., 



(30) r^ = 2 7r5^5 i (q)s J ·(q)(n q + fl(±))í(ü; q -a;) 



i / d 3 qgM)g j (q)(n q +9(±))S(ncq-Tioj) (34) 

Thanks the axial symmetry of the problem and the delta 
function of energy conservation, the three-dimensional in- 
tegral over q in (Q) is better approached in polar coordi- 
nates: d 3 q = q 2 dipd(cos -d)dq. One obtains an expression 
proportional to 

Q 2 t 2 COs(Qtz i: j) 1-t 2 . 2 7TÍ 

dte ,„Z — ^ — ^— sm^f— — 1. (35) 



[(ç/Q) 2 -í 2 ] 2 t 2 



l q/Q' 



with q = qo, a = ao, Q = E /Tic. Moreover, Zij = a(i — j) 
is the distance between i-th and j-th QD's. The crucial 
point is to observe that, for Q/aJ^ 2 = Aj^/A|| , (A|| ~ Q^ 1 ) 
large enough, this integral is dominated by contributions 
around t = cos d = 1; therefore 



r^ = rit J cos[Q % ]. 



(36) 



— 1/2 

Recalling that X± = a ' is the typical length scale of 
carrier confmcmcnt in the x-y planc, this behaviour is 
easily understood: due to the energy-conservation con- 



straints (q^ + q 2 



Q ), for delocalized in-plane 



wavefunctions (with respect to the length scale Ay,) the 
significant fluctuation of q in the considered state is 
small; therefore q z ~ Q. In other words, due to the ex- 
ponential suppression — in the overlap integral — of the 
contributions from phononic modes with non-vanishing 
in-plane components the system behaves as in the pres- 
ence of a single effective phonon mode along the z axis 
resonant with the qubit excitations. As clearly confirmed 
by our numerical analysis reported in Sect. IV, this is an 



extremely important feature of the semiconductor model 
considered: in spite of its 3d nature and of the presence 
of a continuum of decoherence-inducing phonon modes, 
in this regime the carrier subsystem experiences an effec- 
tively 1-d coherent environment, that in a good approx- 
imation can be descxibed by the Circular Model (CM) 
analysed in App. A.o 

This model, parametrized by the dimensionless quan- 
tity Q = Q a, represents a non trivial cxample of a 
register-environment coupling that admits a rich struc- 
ture (as a function of Q) of sub-decoherent encodings. 
From this point of view, it realizes a generalization of 
the replica symmetric model (pure collective decoher- 
ence) discussed in [ï^, that is recovered for Q = 0. Here, 
we limit ourselves to summarize the main result: 
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Safe quantum encoding are possible for the models 
such that é 1 ® is a 4-th roots of the unity, the most effi- 
cient case being the points Q = 0, mod 7r; when all the 
register celis feel the same external coupling the dynam- 
ics is maximally collective thanks the full permutational 
symmetry. 

The existence of infinitely many "màgic" points is 
clearly due to the unphysical nature of the CM that 
allows for undamped interactions between objects sep- 
arated by arbitrary large distances. In realistic systems 
(as the ones investigated in this paper) the cosine de- 
pendence of the T matrix can be only approximated and 
the periodicity with respect to the celi distance eventu- 
ally destroyed by some overimposed decay. In a way, 
the present situation is very similar to having a string of 
(two-level) àtoms in a cavity. coupled with a single reso- 
nant electromagnètic modeEZl. 



IV. SIMULATION OF SUB-DECOHERENT 
DYNAMICS IN A QD ARRAY 

In this section we will present our numerical analy- 
sis of subdecoherent quantum encoding for realistic QD 
structures. 



A. Carrier-Phonon Scattering in a Single QD 
Structure 



As a starting point, let us discuss the role of carrier- 
phonon interaction in a single QD structure. Figure |ï| 
shows the total (emission plus absorption) carrier-phonon 
scattering rate at low temperature (T = 10 K) as a func- 
tion of the energy spacing E for three different vàlues of 
the GaAs QW width (d — 3,4, andònm). Since the en- 
ergy range considered is smaller than the optical-phonon 
energy (36meV in GaAs ), due to energy conservation 
scattering with LO phonons is not allowed. Therefore, 
the only phonon mode A which contributes to the rate 
of Fig. [y is that of acoustic phonons. Again, due to en- 
ergy conservation, the only phonon wavevectors involved 
must satisfy |q| = E/hc s = q, c s being the GaAs sound 
velocity. It follows that by increasing the energy spac- 
ing E the wavevector q is increased, which reduces the 
carrier-phonon coupling entering in the electron- phonon 
interaction and then the scattering rate. This well- 



established behaviour, known as phonon bottleneckxrl is 
typical of a quasi-OD structure. As shown in Fig. for 
E = 5 meV — a Standard value for many state-of-the-art 
QD structures — the carrier-phonon scattering rate is al- 
ready suppressed by almost three orders ofjnagnitude 
compared to the corresponding bulk valuesfl'B 

In addition to the bottleneck scenario discussed so far, 
for a given value of the energy spacing E we see that 
for small vàlues of d we have an increase of the carrier- 
phonon rate. In spite of the reduction of the 3D vol- 



ume available to the carrier states, the overall coupling 
is increased, basically due to the progressive relaxation of 
momentum conservation along the growth (z) direction. 



B. Short-Time analysis 

We will now show that by means of a proper informa- 
tion encoding, i.e., a proper choice of the initial multi- 
system quantum state, and a proper design of our QD 
array, we can strongly suppress phonon-induced decoher- 
ence processes, thus further improving the above single- 
dot scenario. To this end, let us consider a four-QD array, 
which is the simplest noiseless qubit register (see Ap p. A) 



From the short-time expansion discussed in Sect. [IA 
we have numerically evaluated the decoherence rate for 
such QD array choosing as energy splitting E — 5 meV 
and QW width d = 4nm (see Fig. |ï]). As initial state 
we have chosen the singlet IV'Ui.o) [see Eq. ||| defined 
by the dimer partition T>\ = {(1,2), (3,4)}. We stress 
that, when the CM approximation (see App. A) is not 
exactly fulfilled, different singlets have different decoher- 
ence rates. Indeed, the larger is the distance Zij between 
the pair elements in the dimer covering, the greater is 
the deviation from the strictly periòdic behaviour. Thus 
from Eq. ( 17) it follows, for instance, that the singlet 
corresponding to the dimer partition T>2 = {(1, 3) (2, 4)} 
has a greater decoherence rate than \4>t>i,o)\ The deco- 
herence rate obtained from our numerical calculation is 
shown as sòlid line in Fig. |](a) as a function of the inter- 
dot distance a. The uncorrelated-dot decoherence rate 
is also reported as dashed line for comparison. As sug- 
gested by the analysis of the circular model presented in 
App. A, in spite of the 3D nature of the sum over q en- 
tering the calculation of the function T^, [see Eq. ([351)1, 
the decoherence rate exhibits a periòdic behaviour over 
a range comparable to the typical QD length seale. In 
the circular-model approximation (and for T = 0) one 

obtains r-f 1 [|V>Di] — 2Fq ' ) [1 — cos(Qa)], from which it 
follows that for a n — 2mr/Q, (n € n the considered 
state is stable. This effect - which would be natural for 
a 1D phonon system - stems from the exponential sup- 
pression, in the overlap integral, of the contributions of 
phononic modes with non-vanishing in-plane component, 
previously discussed. This 1D behavior is extremely im- 
portant since it allows, by suitable choice of the inter-dot 
distance a, to realize a symmetric regime in which all the 
dots experience the same phonon ficld and therefore de- 
cohere collectively. Figure |2](a) shows that for the partic- 
ular QD structure considered, case C should correspond 
to a decoherence-free evolution of a singlet state, which 
is not the case for A and B (see symbols in the figure). 

In order to better understand how this sort of effective 
1D behaviour depends on the material parameters con- 
sidered, we have repeated the subdecoherence analysis 
of Fig. ^(a) by artificially increasing the GaAs effective 
mass. More specifically, Figs. |2](b) and ||(c) present the 
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same decoherence analysis, respectively, for vàlues of 5 
and lOm*. As we can see, by increasing the efïective-mass 
value the 1D character in (a) is progressively suppressed. 
This can be clearly understood as follows: the increase 
of the effective mass leads to a stronger and stronger 
localization of the 2D harmonic-oscillator wavefunctions 
which, in turn, can easily interact with transverse (xy) 
phonon modes q. 

As far as the unitary component of the Liouvillian is 
concerned, one can easily show that [for any \ip) eigen- 
stateof S z ] F{t) = e~ ltH '^ |V)| 2 = l-(í/r t/ ) 2 +o(í 3 ), 
where 



= 5H\ |V) - (V'l ÓH n W 



(37) 



Figure || shows r^ 1 [IV'üJ a s a function of the inter-dot 
distance a. 

We find an oscillatory behaviour similar to that of 
Fig. ||(a); it stems from the fact that (for the material 
parameters considered) A^- ~ T^ sm[Q(i — j)a + (p], with 
ip <C 7r/2. Thus, for vàlues of a corresponding to a sub- 
decoherent dynamics [see point C in Fig. |](a)] the A 
contribution, also known as polaronic shift, is negligible 
as well. 



C. Time-dependent solution of the Master Equation 

In order to extend the above short-time analysis, we 
have performed a direct num erical integration of the Mas- 
ter equation (see Sect. HA), thus obtaining the reduced 
density matrix p as a function of timc. Also the Lamb- 
shift terms discussed in Sect. || have been taken into 
account. Starting from the same GaAs QD structure 
considered so far, we have simulated the above noise- 
less encoding for a four-QD array. Figure |J shows the 
fidelity as a function of time as obtained from our nu- 
merical solution of the Master equation. In particular, 
we have performed three different simulations — for the 
same initial state |V'i'i,o) — corresponding to the differ- 
ent vàlues of a depicted in Fig. [|(a). Consistently with 
our short-time analysis, for case C we find a strong sup- 
pression of the decoherence rate which extends the sub- 
nanosecond time-scale of the B case (corresponding to the 
uncorrelated-dot rate) to the microsecond time-scale. 

An other quantity which properly describes the 
environment-induced corruption of information is the lin- 
ear entropy 6[p] = tr(p — p 2 ). Its production rate is also 



directly connected to t\ ; indeed for an initial pure prepa- 
ration we have 5(t) = 2 t/r\ +o(í 2 ). for intial pure prepa- 
rations. The time evolution of the linear entropy, as ob- 
tained from our numerical solution of the ME, is reported 
m Fig. | We can clearly recognize an initial transient (of 
the order of T\) in which the register, getting entangled 
with the environment, decoheres; this is followed by a 
subsequent slower relaxation dynamics. 



The time-dependent analysis of Figs. || and || confirms 
that by means of the proposed encoding strategy one 
can realize a decoherence-free evolution over a time-scale 
comparable with typical recombination times in semi- 
conductor materials!! 



V. SUMMARY AND CONCLUSIONS 



We have investigated a possible semiconductor-based 
implcmcntation of the subdccoherent quantum-encoding 
strategy, i. e. error avoiding, recently proposed in |Ï2|. The 
goal is the suppression of phase-breaking processes in a 
quantum register realized by the lowest-energy charge ex- 
citations of a semiconductor QD arrayliia In this case, the 
primary noise source is given by electron-phonon scatter- 
ing, which is considered to be the most efficicnt decoher- 
ence channel in such a systemJ3'EI 

The main result is that, in spite of the 3D nature of 
carrier-phonon interaction in our QD structure, by means 
of a proper quantum encoding as well as of a proper tai- 
loring of the semiconductor structure, one can in princi- 
ple increase the coherence time by several orders of mag- 
nitude with respect to the bulk value. This would allow 
to realize a coherent quantum-mechanical evolution on 
a time-scale longer compared to that of ultrafast optical 
spectroscopy. From this point of view this result might 
constitute an important step toward a solid-state imple- 
mentation of quantum computers. On the other hand, it 
certainly represents a first non-trivial example of a solid- 
state quantum system for which one can apply quantum 
error avoiding strategies. 



As aheady discussed in Sect. III carrier-phonon scat- 
tering is not the only source of decoherence in semi- 
conductors. In conventional bulk materials also carrier- 
carrier interaction is found to play a crucial role. How- 
ever, state-of-the-art QD struc-tures — often referred to 
as semiconductor macroatomsEül — can be regarded as 
few-electron systems basically decoupled from the elec- 
trònic degrees of freedom of the environment. For the 
semiconductor QD array considered, the main source of 
Coulomb-induccd "noise" may arise from the inter-dot 
coupling. However, since such Coulomb coupling van- 
ishes for large vàlues of the QD separation and since the 
proposed encoding scheme can be realized for vàlues of a 
much larger than the typical Coulomb-correlation length 
(see Fig. |^), a proper design of our quantum register may 
rule out such additional decoherence channelsHa 

The actual implementation of the suggested encoding 
relies, of course, on precise quantum state synthesis and 
manipulat ions. This further step, not addressed in this 
paper, represents the most challenging open issue con- 
cerning the ultimate usefulness of the proposed coding 
strategy. 
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FIG. 1. Carrier-phonon scattering rate for a single QD 
structure as a function of the energy splitting E for different 
vàlues of the QW width d at low temperature (see text). 



FIG. 2. (a) Phonon-induced decoherence rate for a 
four-QD array (sòlid line) as a function of the inter-dot dis- 
tance a compared with the corresponding uncorrelated-dot 
rate (dashed line); (b) Same as in (a) but with an artificial 
effective mass of 5m* ; (c) Same as in (a) but with an artificial 
effective mass of lOm* (see text). 



FIG 3. Tjj [IV^i)] as a function of a (see text). 



FIG 4. Fidelity F as a function of time as obtained from 
a direct numerical solution of the Master equation for the 
relevant case of a four-QD array (see text). 



FIG 5. Linear entropy as a function of time as obtained 
from a direct numerical solution of the Master equation for 
the relevant case of a four-QD array (see text). 
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APPENDIX A: CIRCULAR-INTERACTION 
MODEL 

This appendix is devoted to the formal analysis of a 
model with periòdic (environment-induced) interactions 



between register celis. We set T^' = cos[Q (i — j ) 
the resulting model will be referred to as the circular 
model (CM). The dimensionless parameter Q is taken to 
be given by the product of a characteristic wave vector 
(corresponding to an effective one-phonon field) times the 
inter-cell distance. The effective Hamiltonian (13) takes 
the form H = £ Q=± H^ a \Q), with 



H 



(a) _ 



(Al) 



where Sq = Y^,j=i e% a f ( a = ^ z )i are tke present 
Lindblad operators. They fulfill the following commuta- 
tion rclations 



[S* Z , Sq] — ± S. 



±2^ 
0' 



(A2) 



For Q = 0mod2 7r one recovers the global sl(2) àl- 
gebra spanned by the S a, s, to which the Sq's are 
connected by the following unitary transformations 

Uq 



exp(i Q 3 of )• Indeed, 



U aQ S a Ul Q (a 



±) (notice that Uq = U-c 



we have Sq 



In terms 

of these unitary transformations and of the Q = Hamil- 
tonian H = r(-) S+ S- + r(-) S~ S+ the CM model 
(Al) reads 



Hn = 2- 



U vQ H U vQ 



(A3) 



From Eq. (A2) it follows that, for any genèric Q, the two 
terms in the above equation do not commute: the model 
is non trivial, i.e., non integrable. 

Next proposition shows that the analytic structure of 
the CM strongly depends on the input parameter Q, for 
particular Q vàlues it is quite simple and its subdecoher- 
ent coding efficiency is optimal. 

Proposition 1 One has the following integrable points 

i) Q = 0mod27r,^ H {a \2ir) = T^S- a S a , replica 
symmetry. 



ii) Q = 7rmod27r,=> íï^tt) = S~ a {n) S a (n) = 
U^H^(2n)W n . 

Ui) if Q — 7r/2, 3/2 7r mod 2 7r one has r|"^ 2n = 

r( Q )(-l)™ andrlf +2n+1 = 0. The odd- and even- 
site sublattices decouple, and for each sublattice 
case ii) is recovered. 

Notice that for cases i) and ii) 2 Q = mod 2 tt; then 
the ?7 = + and n = — terms in (A3) are identical; the 
model is then unitarily equivalent to the Q — case. The 
latter is clearly diagonalized in the S 2 , S z eigenbasis and 
its spectrum is given by E = X) Q =± ( J, M, r) where 

( J, M, r) = r( Q ) [J ( J + 1) - M (M + a)], (A4) 

J = J mm , . . . , N/2; M = -J, . . . , J; r = 1, . . . , n(J, N), 
in which J m in = ( J m in = 1) for N even (odd), and 
n(J,N) denotes the multiplicity of the sl(2) representa- 
tion labelled by JU 



n(J, N) 



NI (2 J+ 1) 



{n/2 + j + iy.(N/2- jy. 



(A5) 



If N is even and < r( + ) < the lowest eigenvalue is 
Eo = with degeneracy n(0, N), the ground-state man- 
ifold being the singlet sector of the global sl(2). At zero 
temperature one has = 0; therefore all the lowest- 
weight sZ(2)-vectors \ J, —J) are ground states of H. 

Let us consider the N-th roots of the unit (with N 
even) 



Z N = {é Q ^/Q j = 



27Tj 

N ' 



j = 0,...,N-l}. (A6) 



This (multiplicative) group is of course isomorphic to 
the (additive) group Z/NZ = {0, . . . , N - 1}; thus we 
shall use the same notation for both. Here Zn is con- 
sidered a subgroup of Sn- The latter as a natural ac- 
tion on Tíii given by the linear extension p: ®jLi\o~j) i— > 

The operators S^ = S a (Q m ), satisfy to the commuta- 
tion rclations Sg,] = Kff S*^ +m , are the sl(2) 
structure constants. They span the (2jy-graded) Lie àl- 
gebra 

A N = span{5^ / a = z, ±, me Z N } S ef sl(2)i. (A7) 

Let A N the Lindblad operators àlgebra for a genèric 
Q, the following proposition gives a characterization of it 
when Q varies. 

Proposition 2 

i) For a genèric Q (i.e., e 1 ® ^ Zn) one has A N = 
An i whereas for é 1 ® £ Zn one finds 

A% n = span{5| p „, S± (2p+1) /pe Z N } 
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ii) A° N Sí A N S «1(2) 
ui) = s Z(2) e e i sZ(2) , where 

AT/2 

sl{2) e = span{^(-l)V 2 «.} a , 
i=i 

W/2-1 

sZ(2) D = span{ £ (-l) J '<r 2 a i+ i}a (A8) 

3=0 

iv) e lQ ^ G Z* N = Z N - Z 4 => dirnA^ = 3 iV/2 
Proof 

One can check that A N — span{5| p g, Sq^ 2 / p e 

.2} if Q is rationally independent from 2 7r the numbers 
e i2pQ ^ e iQ(2p+i) densely fill the unit circle, from which 

An C -4^Ej. Points ii)-iii) follow from prop. 1, and iv) 
can be checked by a simple calculation □. 

Notice that e l7T / p £ Z N iiï N = 0mod2p( í > = 1,2). 
Remembering that G Ker Hq is annihilated 

by aZZ the generators of — and then that the smaller 
is the àlgebra the greater is the code — Proposition 2 
seems to indicate that the "màgic" Q's possibly relevant 
for subdecoherent encoding are just the ones such that 
e í( 3 G Z N . 

Let us now consider the Q dependence of the symmetry 
structure of our model. 

Lemma Let Gq the (maximal) symmetry group of Hq, 
one has: 

i) Go = Gir — Sn, Ü) Gtt/2 — Gzn/2 = S N / 2 X S N / 2 , Üï) 

for e íQ G S 1 — Z4 one has Gq — Zn- 

Pictorially one can say that in the CM the regis- 
ter has a regular polygon topology that for the special 
points Q — 0, 7r (Q — 7r/2, 3/2tt) collapses to a point 
(dimer) gaining in this way a larger permutational sym- 
metry. This dynamical clustering is asapciated with a 
greater sub-decoherent coding efnciencyO Next proposi- 
tion summarizes in a formal manner the present situa- 
tion. 

Proposition 3 Let N even then 
i) e lQ £Z N ^d N (Q) >0 
ii) d N (Qo) = d N (Q N / 2 ) = n(0, N). 



iü) IfN = 0mod4 => d N (Q N/4 ) = d N/2 (Q ) 2 

iv) e l ®i G Z N =>d N (Qj) = 1. 

v) When e l ® j G Z N the null space is spanned by the 
vectors 

m = ®^i 2 (ioi) - (-iy iio)Wív/2 (A9) 

vi) Let N odd then d^{Q) = OVQ. 
Proof 

The cases e í( 2 = ±1 are isomorphic and have been 
previously discussed. Notice that, if H > one has 
{Í>\H\iP)_ = H\if>) = 0. Moreover, |^) = \f>, Jj) 
where T> is the unique dimer partition of the array with 
\l - k\ = N/2 and 7^, k) = j (mod2) V(Z, k) G V. From 
the second of Eqs. ( [17| ) one finds that 

2 N/2 

(^\H Qj \^) ~ 1 - - ^(-1) J cos(Trj) = 
1=1 

from which the sufficiency parts of i) and v) follow. 
jf g«Q ^ 2, N from prop. 2 and 3, one has that if 
\ip) G Ker Hq then it is in the singlet sector of .4.^ ( prop. 
2) Since C 2 is an irreducible (non-trivial) representation 
space of An such a sector is empty. Points ii)-iii) follow 
directly from prop. 2, 3. Since the S±q's transform ac- 
cording 1-d Ç?Q-irreps, from representation theory it fol- 
lows that dN(Q) (i.e., multiplicity of the 1-d „4^-irrep 
) is equal to the dimension of an irrep of the symmetry 
group Gq - But for e lQ G Z* N one has Gq = Zn, (abeliaal 
therefore its irreps are 1-d, from which point iv) follows.Ell 
Finally, vi) simply stems from the fact that the necessary 
condition S z \ip) — cannot hold for odd N. □ 
To understand in a more constructive fashion, why the 
\ipj) are the (only) subdecoherent states for é l( ^ j G Z* N 
let us consider the following state \ip) G C ( C the 
(global) sZ(2)-singlet sector) such that i) U 2 q \ip) — 
Then U Q \ip) = U Q \tp) = U- Q \ip) = This means 
W) S C\ a= ±U a QC it follows that is annihilated by 

S ±Q = U ±Q Sa U ±Q> ( a = z ' ± )- and therefore by Hq It 
is now easy to check that the states \tpj) of prop. 4 are 
just L r Q j \ f>, 0), the dimer partition V being the only one 
allowing for condition i) to be fulfillcd. 
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P. Zanardi and F. Rossi: Fig. 1 




P. Zanardi and F. Rossi: Fig. 2 
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P. Zanardi and F. Rossi: Fig. 3 




P. Zanardi and F. Rossi: Fig. 4 
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P. Zanardi and F. Rossi: Fig. 5 



